Electronic structures and ferromagnetism in transition metals codoped ZnO 
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We have investigated electronic structures and magnetic properties of potential ZnO based diluted 
magnetic semiconductors: (Fe, Co) and (Fe, Cu) codoped ZnO. The origins of ferromagnetism are 
shown to be different between two. (Fe, Co) codoped ZnO does not have a tendency of Fe-O-Co 
ferromagnetic cluster formation, and so the double exchange mechanism will not be effective. In 
contrast, (Fe, Cu) codoped ZnO has a tendency of the Fe-O-Cu ferromagnetic cluster formation 
with the charge transfer between Fe and Cu, which would lead to the ferromagnetism through the 
double-exchange mechanism. The ferromagnetic and nearly hah- metallic ground state is obtained 
for (Fe, Cu) codoped ZnO. 

PACS numbers: 75.50.Pp, 71.22.+i, 75.50.Dd 
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I. INTRODUCTION 

Electronics utilizing the spin degree of freedom of 
electrons, namely "spintronics" , becomes an important 
emerging field. The spintronics is expected to overcome 
the limits of traditional microelectronics, such as non- 
volatility, data processing speed, electric power consump- 
tion, and integration densities^. Diluted magnetic semi- 
conductors (DMSs) will play a core role in spintronics 
as semiconductors do in electronics due to easy integra- 
tion into existing electronic devices. Two types of DMS 
families have been well studied: II- VI type such as Mn- 
doped CdTe and ZnSe-^ , and III-V type such as Mn-doped 
GaAs^. In particular, the latter attracts great attention, 
because it becomes a ferromagnetic (FM) DMS having 
the Curie temperature Tc ^ llOK. Recent research 
effort has been focused on developing a new FM-DMS 
operating at room temperatur8^i^*SiLSiS. 

Along this line, attempts have been made to fabricate 
ZnO based DMS. ZnO is a wide-gap (Eg ~ 3.44eV) II- 
VI semiconductor, and so it can be used for ultraviolet 
light emitting devices. Jin et almSi fabricated 3d transi- 
tion metal (TM) doped epitaxial ZnO thin films using 
the combinatorial laser molecular-beam epitaxy method. 
However, they have not detected any indication of ferro- 
magnetism. In contrast, Ueda et ali^ observed the FM 
behaviors in some of the Co-doped ZnO films made by us- 
ing the pulsed-laser deposition technique with Tc higher 
than the room temperature. The reproducibility, how- 
ever, was less than 10%. Hence the realization of the FM 
long range order in Co-doped ZnO films is controversial. 

On the other hand, there were also trials to make ZnO 
based DMS by co-doping two TM elements: (Fe,Co) or 
(Fe,Cu). Zni_x(Coo.5Feo.5)xO films were fabricated by 
using the reactive magnetron co-sputtering techniquei^, 
which seemed to have the single phase of the same 
wurtzite structure as pure ZnO up to x=0.15. The room 
temperature ferromagnetism was observed and the rapid 
thermal annealing under vacuum leads to increases in Tc, 
magnetization, and the carrier concentration. Also Cu- 
doped Zni_xFexO bulk samples were fabricatedi^i. The 



bulk sample has the advantage of insensitivity to the de- 
tailed process conditions, over the film sample fabricated 
under the nonequilibrium condition. They observed the 
FM behaviors with Tc~550K and the saturation mag- 
netic moment of 0.75/iB per Fe in Zno.94Feo.05Cuo.01O. 
The saturation magnetic moment increases, as the Cu 
doping ratio increases up to 1%. In addition, the large 
magnetoresistance was observed below lOOK. 

Motivated by these reports of ferromagnetism in TM 
codoped ZnO, we have studied electronic structures and 
magnetic properties of (Fe, Co) and (Fe, Cu) codoped 
ZnO: Zno.875(Feo.5Mo. 5)0.1250 (M=Co or Cu). We have 
used the linearized muffin-tin orbital (LMTO) band 
method in the local spin-density approximation (LSDA). 
To explore the Coulomb correlation and the spin-orbit 
(SO) effect in TM-doped ZnO, we have also employed 
the LSDA+L/(+SO) method incorporating the Coulomb 
correlation interaction U and the SO interactior^. ZnO 
has the wurtzite structure in which anions and cations 
form hexagonal close-packed lattices. The wurtzite ZnO 
is composed of tetrahedrons formed by four O anions. 
For Zno.875(Feo.5Mo. 5)0.1250, we have considered an or- 
thorhombic supercell containing sixteen formula units in 
the primitive unit cell by replacing two Zn atoms by Fe 
and Matoms (Zni4FeiMiOi6). For the lattice constants, 
we assumed those of pure ZnO with a — 6.4998, b = 
11.2580, c= 5.2066 A. 



II. ZnO, Fe- and Co-doped ZnO 

First, we have checked the electronic structure of pure 
wurtzite ZnO without doping elements. The overall band 
structure of the present LMTO result is consistent with 
existing resultsi^iiSiii. As usual in the LSDA calcula- 
tions, the obtained energy gap ~ 0.7 eV is only about 
20% of the experimental value. Also the position of Zn 
3d band (—3.0 ^ —5.0 eV with respect to the valence 
band top) is much shallower than the Zn 3d spectrum ob- 
tained by photoemission experiment (~ —9.0 eV)^^. The 
LSDA-f [/ band method improves the results of LSDA, 
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FIG. 1: The LSDA+C/ band structure of ZnO with (7=3. OeV 
for Zn 3d electrons. The lowest and intermediate bands cor- 
respond to mainly Zn 3d and O 2p bands, respectively. Zn 3d 
bands in the LSDA are too shallow to become mixed with O 
2p bands. 



but not SO satisfactorily. As shown in Fig. ^ the energy 
gap is increased to ~ 1.0 eV, and the position of Zn 'id 

band becomes deeper (-6.0 7.0 eV). The LSDA+C/ 

result for Zn 3d position is consistent with result obtained 
by the GW band calculationiiS*ii. Still the energy gap 
and Zn id position are smaller and shallower than exper- 
imental values. 

We then have performed the LSDA band calculations 
for single TM-element doped ZnO: Fe- and Co-doped 
ZnO with 6.25% concentration of TM element. Due to 
heavy computational load, we have not considered the 
Coulomb correlation U for Zn 3d electrons in the super- 
cell calculations of doped ZnO systemsi^. Note that elec- 
tronic structures of TM doped ZnO have already been 
reported for 25% concentration of TM element by us- 
ing the Korringa-Kohn-Rostoker band method combined 
with the coherent potential approximation^. They found 
that V-, Cr-, Fe-, Co-, and Ni-doped ZnO would have 
the FM ground states rather than the spin- glass states. 
With this background, we thus consider below only the 
FM states for Fe-, and Co-doped ZnO with more realistic 
TM doping concentration. 

Figure|21shows the LSDA total density of states (DOS) 
and TM 3d projected local DOS (PLDOS) for Zni^^Fe^O 
and those for Zni_xCoxO (x — 0.0625). For both Fe- and 
Co-doped ZnO, we have obtained nearly half -metallic 
electronic structures, that is, the conduction electrons at 
the Fermi level Ep are almost 100% spin polarized. The 
minority spin 3d states of both Fe and Co near Ep are 
seen to be hybridized slightly with the conduction band 
(see Fig.O. Both for Fe-doped and Co-doped ZnO, the 
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FIG. 2: The LSDA total and TM 3d PLDOSs of 
Zno.9375Feo.o625 and Zno.g375Coo.o625 0. 



Fermi level cuts the sharp minority spin eg states. Note 
that Fe and Co are located at tetrahedral centers formed 
by O ions, and so eg states are lower in energy than 
t2g states. Since there is nearly one more d-electron in 
Co-doped ZnO, the Fermi level is located near the valley 
between Co eg and t2g minority spin state. The exchange 
splittings are larger than the crystal field splittings, re- 
flecting high spin states of Fe and Co in ZnO. Especially, 
the exchange splitting in Fe-doped ZnO is very large to 
locate both eg and t2g states deep in energy, and so the 
majority spin Fe t2g states become fully hybridized with 
O 2p states to yield a broad band. In contrast, the t2g 
states in Co-doped ZnO are shallow and located above 
the O 2p valence band and so the hybridization becomes 
weak. 

The total magnetic moments 4.44 and 3.25 /xb for 
Fe- and Co-doped ZnO come mostly from Fe (4.11 ^b) 
and Co (2.92 ^b) ions, respectively. These results sug- 
gest electron occupancies of d^ (Fe^+) and d'' (Co^+), 
respectively. Figure 13 provides the band structure of 
Zno.9375Coo.o6250 near Ep. Note that the size of circle 
represents the amount of Co 3d component in the wave 
function. It is seen that the rather flat majority spin Co 
3d states are located by ~ 1.0 eV below the conduction 
band consisting of mainly Zn 4s states. The minority 
spin Co 3d states are located near and above Ep man- 
ifesting hybridization with the Zn 4s conduction band 
states. However, we do not expect that either Fe- or Co- 
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FIG. 3: Band structure of Zno.9375Coo.o625 near the Fermi 
level (the upper panel for the majority and the lower for the 
minority spin band). The size of circle represents the amount 
of Co 3d component in the wave function. 



doped ZnO in nature has a stable metallic ferromagnetic 
ground state, as obtained above. The high DOSs at Ep 
for both Fe- and Co-doped ZnO would drive the possible 
structural instability or become reduced substantially by 
the Coulomb correlation interaction between TM 3d elec- 
trons. Indeed, LSDA-|-[/ band calculation with U =5.0eV 
for Co 3d electrons yields the insulating ground state for 
Co-doped ZnO, distinctly from the LSDA band resultaSS. 



III. (Fe,Co) codoped ZnO 

Now we have performed the LSDA band calculations 
for (Fe, Co) codoped ZnO: Zno.875(Feo. 5000.5)0.1250. We 
have examined magnetic properties of supercell by vary- 
ing the separation between Fe and Co: 3.2499, 5.6055, 
and 6.4998 A. 3.2499 A corresponds to a nearest Fe-Co 
separation with Fe-O-Co configuration in the a-b plane 
of the orthorhombic supercell, while 5.6055 and 6.4998 A 
to farther Fe-Co separations with Fe-O-Zn-O-Co configu- 
rations along the diagonal direction and in the a-b plane, 
respectively. Total energies are nearly the same among 
three cases: the shortest 3.2499 A case has the lowest 
total energy by only ~ 3 mRy. This result reflects that 
there will not be any noticeable TM clustering effect in 
the (Fe, Co) codoped ZnO. Further we have found that 
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FIG. 4: The LSDA total and TM 3d PLDOSs of 

Zno. STsFeo. 0625 Coo. 0625 O. 



the FM configuration of Fe and Co spins are found to be 
slightly more stable than the antiferromagnetic (AFM) 
configuration for all three cases. 

Figure H shows the LSDA DOS of Zni^axFe^Co^O 
(x = 0.0625) for a Fe-Co separation of 5.6055 A. The 
3d PLDOS in the codoped case looks as if it is just a 
simple sum of Fe and Co PLDOSs of Fig. [3 There 
is no indication of charge transfer between Fe and Co, 
and so Fe^+ and Co^+ valence configurations are re- 
tained in the (Fe, Co) codoped ZnO. Hence most prop- 
erties are similar to each Fe- and Co-doped case, such 
as the degree of hybridization, nearly half-metallic na- 
ture, and the local magnetic moments. This suggests that 
the double-exchange mechanism^^ will not be effective in 
Zni_2xFexCoxO, because the kinetic-energy gain through 
the hopping of spin-polarized carriers between Fe and Co 
ions does not seem to occur. Thus to explain the observed 
FM ground state in Zni_2xFexCoxO, one needs to invoke 
another exchange mechanism between Fe and Co, such 
as the RKKY-type exchange interaction mediated by Zn 
4s carriers or conduction carriers induced by oxygen va- 
cancies. However, one cannot rule out the formation of 
separated Fe or Co metallic clusters in Zni_2xFexCoxO 
which would exhibit the ferromagnetism. Also the possi- 
ble formation of impurity phases such as spinel CoFe204 
is worthwhile to be checked carefully. These features re- 
main to be resolved more in the experiments. 
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FIG. 5: The LSDA and LSDA+(/+SO {U=3.0 eV) total and 
Cu 3d PLDOSs for Zno.9375Cuo.o625 0. 



IV. Cu-doped ZnO 

Before discussing (Fe, Cu) codoped ZnO, we have 
examined electronic structure of Cu only doped ZnO: 
Zno.9375Cuo.o6250. Interestingly, as shown in Fig. we 
have obtained the stable FM and half-metallic ground 
state for Cu-doped ZnO in the LSDA. The total mag- 
netic moment is 1 /is and the local moment of Cu is 0.81 
Hb, corresponding to the Cu^+ (cP) valence state. The 
Cu 3d PLDOS in Fig. EJb) indicates that Cu has the 
DOS of an intermediate spin state. In contrast to Fe- 
or Co-doped case, the empty Cu 3d states are located in 
the gap region without the hybridization with the Zn 4s 
conduction band, and thus Cu 3d states are strongly lo- 
calized. We presume that this would be the reason why 
the Cu solubility in ZnO is so low: ~ 1% at most. 

Since the Cu 3d states near Ep correspond to partially 
occupied atomic-like t2g states, Cu ions would have large 
orbital magnetic moment. Indeed the LSDA-|-t/-|-SO cal- 
culation yields the substantial orbital magnetic moment 
of l.ObfiB with the insulating electronic structure^S. The 
large orbital moment arises from occupied minority spin 
t2g states split by the Coulomb correlation and the spin- 
orbit interaction^^. The orbital moment is polarized in 
parallel with the spin moment, and so the total magnetic 
moment amounts to 1.85/^s/Cu. 

In general, an ion at the tetrahedral center with low 
spin cP state would be Jahn- Teller active. In fact, there 
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FIG. 6: The LSDA total and TM 3d PLDOSs for 

Zno.STsFeo. 0625 Cuo. 0625 O. 



was a report that the account of the dynamical Jahn- 
Teller effect is necessary to explain the observed param- 
agnetic susceptibility for Cu-doped Zn024. To examine 
the Jahn- Teller effect in Cu-doped ZnO, we have con- 
sidered the local distortion of a tetrahedron around Cu 
by tilting oxygen ions at the corners. The oxygen ions 
are tilted by 13° toward the xy-plane with retaining 
Cu-0 separations (here xyz coordinates represent the lo- 
cal principle axis). The LSDA band calculation for this 
system yields a more stable insulating ground state with 
enhanced spin magnetic moment of 0.99/ib/Cu. That is, 
due to the Jahn- Teller distortion, Cu 3d electrons become 
more localized and the system becomes insulating. The 
orbital moment in this case would be quenched due to 
the Jahn- Teller effect. 



V. (Fe,Cu) codoped ZnO 

As is done for (Fe, Co) codoped case, we have studied 
magnetic properties of (Fe, Cu) codoped Zni_2xFexCux0 
(x = 0.0625) with varying the separation between Fe and 
Cu: 3.2499, 5.6055, and 6.4998 A. For all three cases, the 
FM configurations of Fe and Cu spins are found to be 
more stable than the AFM configurations, as for the (Fe, 
Co) codoped case. In contrast to the (Fe, Co) codoped 
case, however, we have found that the shortest Fe-Cu 
configuration becomes much more stable with respect to 
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other two cases by ~ 30 mRy. This total energy differ- 
ence is very large when comparing with ~ 3 mRy differ- 
ence for the (Fe, Co) codoped case. This result indicates 
that Fe and Cu ions in (Fe, Cu) codoped ZnO have a 
tendency to form the Fe-O-Cu clusters. For 1 % Cu- 
doped Zno.95Feo.05O bulk sample, Han et al}i observed 
rather small saturated magnetic moment of 0.75 fiB per 
Fe. This value corresponds to only about 1/5 of ideal 
Fe^"*" local moment A hb- Incidentally, the value 1/5 is 
matched with Cu/Fe doping ratio, suggesting a possibil- 
ity that the only Fe's forming the Fe-O-Cu clusters would 
give rise to the FM moment, that is, other Fe ions do 
not produce the long range FM order. As shown below, 
the FM ground state in (Fe, Cu) codoped ZnO can be 
understood based on the enhanced double-exchange-like 
interaction through the Fe-O-Cu clustering effect. 

Figure El presents the LSDA DOS of 
Zno.875Feo.o625Cuo.o6250 for the shortest 3.2499 A 
configuration. Above all, one can notice the strong 
hybridization between Fe and Cu states. Further, it 
is clearly seen in Fig. |SJc) that some new states emerge 
in the Cu 3(i minority spin PLDOS between eg and 
t2g. The new states have dxy-\\\ie characters directing 
toward Fe ionsS^. Cu "id PLDOS has a reduced exchange 
splitting, manifesting the DOS characteristic of the 
low spin state. This feature is different from other two 
configurations with larger Fe-Cu separation, which have 
an intermediate spin state as in Cu only doped ZnO. 
For all three cases, the hybridization between Cu 'id 
and the conduction band exists, which is again distinct 
from Cu only doped ZnO. Due to this hybridization, 
the conduction carriers of mostly Zn 4s states become 
reduced. Noteworthy from PLDOSs of Fig. |H1 is that 
there occurs charge transfer from Fe to Cu, and accord- 
ingly Fe and Cu are likely to have nominal Fe'^"'" (d^) 
and Cu^^ {d^'^) configurations, respectively. As a result, 
the electron occupancy at Cu site increases, and so Cu 
has the reduced spin magnetic moment of 0.51 as 
compared to 0.81 ^b in Cu only doped ZnO. This feature 
explains the experimental result of the reduced number 
of carriers in (Fe, Cu) codoped ZnO with respect to that 
in Fe only doped ZnOi^. The charge transfer from Fe 



to Cu is expected to disturb the Jahn- Teller distortion 
at Cu sites and concomitantly make a system metallic. 
Further, it will cause the mixed-valent occupancies for 
Fe (Fe2+-Fe3+) and Cu (Cu2+-Cui+) ions, and the 
consequent double-exchangc-likc interaction is expected 
to induce the ferromagnetism in (Fe, Cu) codoped ZnO. 



VI. CONCLUSION 

We have investigated electronic structures of (Fe, Co) 
and (Fe, Cu) codoped ZnO together with those of Fe- 
, Co-, and Cu-doped ZnO. We have also explored ori- 
gins of observed ferromagnetism in (Fe, Co) and (Fe, Cu) 
codoped ZnO. The single TM-doped systems would not 
have stable metallic ferromagnetic ground states, due to 
large Coulomb correlation or other structural instability 
effects in Fe- and Co-doped ZnO and the Jahn- Teller ef- 
fect in Cu-doped ZnO, respectively. For (Fe, Co) codoped 
ZnO, we have found no indication of charge transfer be- 
tween Fe an Co, suggesting that the double-exchange 
mechanism will not be effective for the observed ferro- 
magnetism in (Fe, Co) codoped ZnO. Therefore one needs 
to invoke another exchange mechanism between Fe and 
Co, or the possible formation of impurity phases is to be 
checked carefully. In contrast, for (Fc, Cu) codoped ZnO, 
there is a tendency to form the Fe-O-Cu clusters and so 
to give rise to charge transfer from Fe to Cu ion. This 
tendency causes the mixed-valent occupancies between 
Fe and Cu, and accordingly the double-exchange-like in- 
teraction is expected to induce the ferromagnetism in 
(Fe, Cu) codoped ZnO. The FM and nearly half-metallic 
ground state is obtained for (Fe, Cu) codoped ZnO. 
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